Local moment fluctuations in an optimally-doped high T c superconductor 
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We present results of neutron scattering experiments on YBa2Cu3 06.95 (T C =93K). Our results 
indicate that magnetic collective modes due to correlated local moments are present both above 
and below T c in optimally doped YBCO. The magnon-like modes are robust and not overdamped 
by itinerant particle-hole excitations, which may point at a substantial static or slowly fluctuating 
charge inhomogeneity. We compare the experimental results to predictions of the Fermi liquid (FL) 
theory in the Random Phase Approximation (RPA). 

PACS numbers: 74.72.-h, 74.25. Gz, 74.20.Mn 



High temperature superconductivity occurs in layered 
copper oxides at compositions between undoped insu- 
lating and strongly doped conventional metallic phases. 
Some propose that exotic states inherited from the an- 
tiferromagnetic parent Mott insulators are essential for 
superconductivity^ whereas others claim that supercon- 
ductors with the highest transition temperatures, T c , 
conform to the standard model of metals, the Fermi liq- 
uid theory^. 

Magnetic spectra of the copper oxides change dramat- 
ically as a function of doping. At x = 0, these com- 
pounds are Mott insulators with strong on-site repul- 
sion that localizes electrons whose spins order antifer- 
romagnetically. Their low energy excitations are spin 
waves (magnons), which disperse upward from the anti- 
ferromagnetic (AF) ordering vector— 7 . Doping disrupts 
long range AF order, but dispersive excitations simi- 
lar to magnons remain in undcrdoped superconductors 
(0 < x < .13) at high energies^i 7 ^. However, such 
magnetic excitations have not so far been detected in 
superconductors with T c > 80K. Previous experimental 
studies of the compounds with the highest T c s focused 
primarily on the resonance feature around 40me V 9 ' 10 i 11 
appearing only below T^. Both Fermi liquid and non- 
Fermi liquid scenarios have been proposed to explain the 
resonance 1 ^^^^. 

Here, we focus on magnetic excitations in optimally 
doped YBa2Cu3 06.95 (T c = 93 K) above T c . Since previ- 
ous studies failed to clearly distinguis h 12 : 18 the magnetic 
signal from the nuclear background, we used much longer 
counting times (more than lhr/point in some scans) to 
significantly reduce statistical error. We report that the 
collective magnetic excitations do exist in YBa2Cu30g.95 
and that their intensity is appreciable in comparison to 
the magnon intensity in the parent insulator and the un- 
derdoped YBa2Cu30-6.6- 

The experiments were performed on the IT triple-axis 



spectrometer at the ORPHEE reactor at Saclay utilizing 
a double focusing Culll monochromator and a pyrolytic 
graphite (PG002) analyzer. The high quality sample of 
YBa2Cu3 6.95 was the same as in an earlier study-i 9 -. 

Following Ref. [2(], magnetic collective modes were sep- 
arated from the nuclear background by subtracting spec- 
tra recorded at a high temperature where the magnetic 
modes are largely suppressed. Fig. QJa) and (b) give ex- 
amples of the raw data at 52.5 and 60 meV where there 
is a clear double peak structure around /i=0.5, the AF 
wavevector of the undoped parent compounds, at 10K 
and 100K. The intensity of this feature is reduced on 
heating above 300K, as expected from magnetic collec- 
tive modes. Some or all intensity that remains at the high 
temperatures may be magnetic but to be on the conserva- 
tive side, we subtract the entire high T spectrum divided 
by the Bose factor from the 10K and 100K data. 

The measured energies were chosen to minimize sys- 
tematic errors due to contamination by nuclear scat- 
tering. The biggest error in our experiment is an un- 
derestimate of the magnetic signal because some of it 
(10-20%) may remain at high temperatures where the 
background was measured. Any broad continuum is in- 
distinguishable from the background and is also not in- 
cluded in the extracted intensities. Fig. [5] shows that 
normal state magnetic intensity extracted by this pro- 
cedure obeys crystal symmetry appearing with a simi- 
lar intensity in the vicinity of equivalent AF wavevec- 
tors Q = (1.5, 0.5, 1.7) and (1.5, -0.5, 1.7). The observed 
magnetic signal was centered at equivalent antiferromag- 
netic wavevectors Q = (1.5,0.5,1.7), Q = (0.5,0.5,5.4), 
and Q = (1.5,1.5,1.7) in units of (27r/a, 2n/b, 2tt/c), 
(where a, 6, and c are lattice constants (data not shown 
for the latter)). Its intensities in different Brillouin zones 
agree with the magnetic form factor confirming the va- 
lidity of this procedure. Since YBCO is a bi-layer com- 
pound, it has magnetic excitations of acoustic and optic 
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FIG. 1: (Color online) Raw data and background subtraction 
procedure. (a,b) Scans at 52.5 and 60meV at 10K, 100K, 
and 340K offset by a constant for clarity. Lines are guides 
to for the eye. Note that the peak in (b) at 340K is most 
certainly nonmagnetic, (c) Scans at 50meV through the anti- 
ferromagnetic wave vector divided by the Bose factor. Solid 
lines represent smoothed data. The feature centered at h=0.5 
at 10K and 100K is magnetic because it is suppressed at 330K 
as the underlying correlations become weaker, leaving behind 
a "hump" of nuclear scattering peaked at h=0A. The back- 
ground levels near h=0 and h=l do not exactly match due 
to small nuclear contributions whose temperature-dependence 
does not follow the Bose factor (e.g. multi-phonon or incoher- 
ent elastic nuclear scattering). Linear corrections were added 
to the 100K and 330K data to match the backgrounds near 
h=0 and h=l with the result shown in (b). These corrections 
were small is all cases (the 50 meV data are the worst case); 
(d) 50 meV scans after adding the linear terms to the 100K 
and 330K data. We assign the intensity difference between 
10K/100K data and 330K curve to magnetic scattering. Er- 
ror bars represent s.d. 



FIG. 2: (Color online) Magnetic intensity at 120K at 55meV 
extracted by the procedure described in the text and Fig. 
Q] The data in the vicinity of equivalent AF wave vectors 
Q=(1.5, 0.5, 1.7) (left panel) and Q=(l. 5,-0.5,1.7) (right 
panel) are shown in arbitrary units but on the same scale. 
Small differences between the two curves are consistent with 
different resolution functions at the two wavevectors and/or 
statistical error. This figure demonstrates that the intensity 
that we assign to magnetic scattering has the correct symme- 
try. 
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character. Our choice of L=1.7 and 5.4 selects acoustic 
magnetic excitations; a separate study will be necessary 
to investigate optic modes. 

Fig. [3] shows the magnon scattering intensity of in- 
sulating YBa2Cu306.i and of the underdoped Tc=60K 
superconductor, YBa2Cu306.6, a t 52.5meV plotted to- 
gether with that of the optimally-doped sample at 100K 
on the same intensity scale. IN the underdoped sam- 
ples it was obtained by subtracting a constant back- 
ground since significant magnetic intensity is still present 
at 300K. To correct for sample volume we scaled to- 
gether spectra of the bond-buckling phonon at E = 42.5 
meV. Its eigenvector is independent of doping because 
it does not contain any chain oxygen vibrations^. The 
phonon was measured in identical spectrometer configu- 
rations (at Q= (0.5, 0.5, 11)) in the two samples (Fig. 
|3Ji). Therefore, no resolution corrections were needed. 
The magnon peak intensity in YBa2Cu30g.i is six times 



FIG. 3: (Color online) Comparison of magnetic scattering 
of the YBCOe.i YBC0 6 . 6 , and YBCOe.gs samples.(a) Spec- 
tra of the Cu-0 buckling mode scaled by the factor equal to 
the ratio of the volumes of the three samples, (b) Magnetic 
scattering intensities plotted on the same intensity scale after 
normalization by the volume ratio extracted from comparison 
of phonon intensities in (a). 



stronger than that of the normal state magnetic signal in 
YBa2Cu30e.95 whereas the peak intensity of the T C =60K 
sample measured at 100K was three times stronger. How- 
ever, the magnetic signal in YBa2Cu30g.95 is signifi- 
cantly broader than in the insulating phase and in the 
underdoped sample (In the insulator it is nearly Q- 
resolution- limited) . Thus their q-integrated intensities 
of the three samples must be of the same order of mag- 
nitude around 50meV at 100K after two-dimensional Q- 
integration. In making this statement, we are assum- 
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FIG. 4: (Color online) Magnetic intensities at different ener- 
gies, (a) Magnon peak in YBa2Cu30e.i at 52.5 meV after 
subtraction of a flat background. Solid line represents the 
spin wave model in Ref. e convoluted with the spectrometer 
resolution, (b-h) Magnetic scattering in YBa2Cu306.95 as a 
function of the distance from the antiferromagnetic ordering 
wave vector. All spectra have been normalized to the same 
scale. A constant background of 530 and 630 counts was sub- 
tracted in (a) and (b) respectively, whereas the other magnetic 
spectra were obtained as described in the text and Fig[T] The 
weak broad 100K signal at 35meV may be an artifact of the 
background subtraction procedure. Error bars represent s.d. 
Our estimate of systemetic errors given by vertical black lines 
is based on the temperature variation of the background. 



cantly less steep than the magnon dispersion in insulat- 
ing YBa 2 Cu306.i. (Fig. 5a, b) Assuming linear extrapo- 
lation, it crosses the zone boundary around 130/100 meV 
in the [100]/ [110] direction, respectively, compared with 
220meV for the magnons in the insulator—. It also dis- 
perses significantly less steeply than the upper branch in 
YBa2Cu306.6: which has a much narrower lineshape at 
52.5meV (Fig. 3) with the magnetic intensity remaining 
near to q=(0.5,0.5)£ 
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FIG. 5: (Color online) Magnetic susceptibility of 
YBa2Cu30e.95. ( a /b) Experimentally measured mag- 
netic susceptibility difference between 10K/100K respectively 
and 300K in the 110 direction after removing effects of the 
resolution. Different ad. hoc. functional forms were tried 
until their convolution with the spectrometer resolution 
agreed with the data, (c) Calculations based on FL/RPA 
(Ref.—) for the difference between the superconducting state 
(10K) and normal state (300K). (d) RPA result for the 
normal state. Here the difference between 100K and 300K 
is negligibly small (in contrast with experiment), so only 
100K result is shown without taking the difference. Note 
the logarithmic intensity scale when comparing experimental 
and calculated spectra. 



ing that magnetic scattering is broad in q and has a 
fourfold symmetry around the reduced antiferromagnetic 
wave vector Qaf— (0.5,0.5) in a twinned sample like ours. 

Fig. d](b-h) shows magnetic spectra of YBa 2 Cu3 06.95 
measured at different energies at 10K and 100K, respec- 
tively, and Fig. [5]^a,b) shows the cut along [110] af- 
ter removing effects of the resolution. It is based on 
the data shown in Fig. H and Ref. [H At 100K we 
observe a branch above 43 meV dispersing away from 
Qaf- Its q-width is 0.3A -1 full width at half maximum 
(FWHM), which corresponds to a correlation length of 
about 20A Upon cooling to 10K, the scattering inten- 
sity below 60 meV increases and a downward-dispersing 
branch with a constant q-integrated intensity appears 
between 33 and 41 meV below T c as shown in Refs. 
[T9II20I . The q-width of the lower branch is relatively nar- 
row (Fig. 5(a)), corresponding to coherent domains of 
about 55A^. The upward-dispersing branch is signifi- 



Fig. [5tc, d) shows the prediction of the FL theory in the 
standard RPA approximation (FL/RPA )13 1 14 1 15 1 16 us i n g 
the calculation based on Refi 1 - 6 -. The RPA expression for 
the spin susceptibility is: 

Xrpa(<1^) = - — (1) 

1 - 3qXo (q> w ) 

where Xq (q, oj) is the bare transverse susceptibility and 
<7 q is an electron-electron interaction (four-point vertex) , 
which in general can be momentum dependen t 13 ' 14 ! 15 ' 16 . 
Once the parameters of the model were picked to repro- 
duce the known Fermi surface, band width of optimally- 
doped copper oxides, the superconducting gap^i, and the 
resonance peak energy, there was no further possibility to 
adjust parameters to alter the calculation results in any 
significant way. Imaginary part of Xrpa^Qt 1 ^) ^ s P ro " 
portional to the neutron scattering cross section and can 
be directly compared with experiment. Below T c sharp 
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downward-dispersing resonance collective modes appear 
in the calculation due to an excitonic effect inside the 
d x 2_ y 2 -wave superconducting gap (Fig. 5c) 16 . This 
feature appears to be in good agreement with the ex- 
periment (Fig. 5 a,c). The calculation also predicts an 
upward-dispersing feature below T c . However it is an 
order of magnitude weaker (relative to the downward- 
dispersing branch) and has a steeper effective dispersion 
then the measured magnetic signal. 

Above T c FL/RPA predicts only a broad temperature- 
independent particle-hole continuum and no collective 
modes (Fig. 5d). Within a simple Fermi-liquid picture 
one expects weak temperature dependence of the Im% 
proportional to the UbT/Ef with Ep being the Fermi 
energy. In optimally-doped YBCO this ratio is too small 
to be detected at the energies that we probed. In con- 
trast with the FL/RPA result, our experiment shows that 
the spin signal grows significantly on cooling from 340K 
to 100K. Including the renormalization of the quasipar- 
ticles due to scattering by the spin fluctuations within 
a so-called fluctuation exchange approximation (FLEX) 
yields a much stronger temperature dependence of the 
spin excitations at low energies but not at the high en- 
ergies that we investigated 22 . It is not entirely clear how 
further renormalization of the interaction strength (ver- 
tex corrections) will affect the results of the FLEX cal- 
culation. Our preliminary estimate is that it will further 
reduce the predicted temperature dependence. 

The above analysis leads us to conclude that the 
FL/RPA picture does explain magnetic intensity above 
~45meV in neither the normal state nor in the super- 
conducting state. The 52.5 and 60 meV data (Figs. 1, 
3c, d) show incommensurate collective modes both above 
and below T c whose dispersion is similar to magnons in 
the parent insulator. (Fig. 5b) (Note that the scans at 
55meV in Fig. 2 do not have a double-peak structure be- 



cause they were centered not at q=(1.5 0.5 1.7), but at 
one of the incommensurate satellites: Q=(1.6 0.5 1.7).) 
Based on this observation we assign these modes to col- 
lective excitations of local moments inherited from the 
insulating phas o 23 i 24 ' 25 i 26 as opposed to itinerant quasi- 
particlcs inherited from the overdoped conventional FL 
phase. These local moments must coexist with the itin- 
erant particle-hole continuum spanned by the Fermi sur- 
face. 

Our results suggest a qualitative similarity between the 
YBCO family and the La2_ a: Sr a ;Cu04 family of high T c 
superconductors. Wakimoto et al>21 found that low en- 
ergy magnetic excitations peaked near Qaf disappear 
upon increased doping whereas the high energy fluc- 
tuations persist far into the overdoped regime 2 ^. In 
YBa2Cu306+x the lowest energy for detecting magnetic 
excitations due to local moments also increases with in- 
creasing doping reaching approximately 45 meV at opti- 
mal doping as shown in our investigation. This apparent 
gap in the measurable magnetic signal probably appears 
because the magnetic modes are not magnons in the strict 
sense but rather precursor fluctuations to the formation 
of antiferromagnetic clusters. Developing a detailed the- 
ory of this phenomenon is outside the scope of our study. 

One possibility that can be ruled out is that our sample 
may still be in the pseudogap state at 100K and thus the 
observed local magnetism is associated with the pseudo- 
gap. Figure 2 shows that the magnetic signal persist at 
least up to 120K, which is above the pseudogap temper- 
ature reported for T C =93K YBCO. Thus the magnetic 
signal above T c at optimal doping is not associated with 
the pseudogap 32 

Acknowledgement: D.R. would like to thank S.A. 
Kivelson, J. Zaanen, and M. Vojta for valuable comments 
on earlier versions of the manuscript. I.E. would like to 
thank T. Dahm for helpful discussions. 



1 P.A. Lee, N. Nagaosa, and X.-G. Wen, Rev. Mod. Phys. 
78, 17 (2006). 

2 Ar. Abanov, A.V. Chubukov, and J. Schmalian, Adv. 
Phys. 52, 119 (2003). 

8 J.M. Tranquada, G. Shirane, B. Keimer, S. Shamoto, and 
M. Sato, Phys. Rev. B 40 4503 (1989). 

4 S.M. Hayden, G. Aeppli, T.G. Perring, H.A. Mook, and F. 
Dogan, Phys. Rev. B 54, R6905 (1996). 

5 M. Arai, T. Nishijima, Y. Endoh, T. Egami, S. Tajima, K. 
Tomimoto, Y. Shiohara, M. Takahashi, A. Garrett, and S. 
M. Bennington, Phys. Rev. Lett. 83, 608 (1999). 

6 J. M. Tranquada, H. Woo, T. G. Perring, H. Goka, G. D. 
Gu, G. Xu, M. Fujita, and K. Yamada, Nature 429, 534 
(2004). 

7 S. Hayden, H.A. Mook, Pengcheng Dai, T.G. Perring, and 
F. Dogan, Nature 429, 531 (2004). 

8 V. Hinkov, P. Bourges, S. Pailhes, Y. Sidis, A. Ivanov, C. 
D. Frost, T. G. Perring, C. T. Lin, D. P. Chen, and B. 
Keimer, Nature Physics 3, 780 (2007). 

9 J. Rossat-Mignod, L. P. Regnault, C. Vettier, P. Bourges, 



P. Burlet, J. Bossy, J. Y. Henry, and G. Lapertot, Physica 
C 185-189, 86 (1991). 

10 H. He, P. Bourges, Y. Sidis, C. Ulrich, L. P. Regnault, 
S. Pailhs, N. S. Berzigiarova, N. N. Kolesnikov, and B. 
Keimer, Science 295, 1045 (2002). 

11 H. F. Fong, P. Bourges, Y. Sidis, L. P. Regnault, A. Ivanov, 
G. D. Gu, N. Koshizuka, and B. Keimer, Nature 398, 588 
(1999); 

12 H.F. Fong, B. Keimer, P.W. Anderson, D. Reznik, F. Do- 
gan, and LA. Aksay, Phys. Rev. Lett. 75, 316 (1995). 

13 E. Dernier, and S.C. Zhang, Phys. Rev. Lett. 75, 4126 
(1995). 

14 D.Z. Liu, Y. Zha, and K. Levin, Phys. Rev. Lett. 75, 4130 
(1995). 

15 F. Onufrieva, and P. Pfeuty, Phys. Rev. B 65, 054515 
(2002). 

16 I. Eremin, D.K. Morr, A.V. Chubukov, K.-H. Bennemann, 
and MR. Norman, Phys. Rev. Lett. 94, 147001 (2005). 

17 M. Vojta, T. Vojta, and R.K. Kaul, Phys. Rev. Lett. 97, 
097001 (2006). 



■5 



Hyungje Woo, Pengcheng Dai, S. M. Hayden, H. A. Mook, 
T. Dahm, D. J. Scalapino, T. G. Perring, and F. Dogan, 
Nat. Physics 2, 604 (2006). 

D. Reznik, P. Bourges, L. Pintschovius, Y. Endoh, Y. Sidis, 
T. Masui, and S. Tajima, Phys. Rev. Lett. 93, 207003 
(2004). 

P. Bourges, Y. Sidis, H. F. Fong, L. P. Regnault, J. Bossy, 
A. Ivanov, B. Keimer, Science 288, 1234 (2000). 
A. Damascelli, Z. Hussain, and Z.-X. Shen, Rev. Mod. 
Phys. 75, 473 (2003). 

T. Dahm and L. Tewordt, Phys. Rev. B 52, 1297 (1995). 
CD. Batista, G. Ortiz, and A.V. Balatsky, Phys. Rev. B 
64, 172508 (2001). 

M.V. Eremin, AA. Aleev, and I.M. Eremin, Zhur. Eksp. 



Teor. Fiz. 133, 862 (2008). [JETP 106, 752 (2008)]. 

25 I. Sega, P. Prelovsek, and J. Bonca, Phys. Rev. B 68, 
054524 (2003). 

26 A. Sherman and M. Schreiber, Phys. Rev. B 68, 094519 
(2003). 

27 S. Wakimoto, K. Yamada, J. M. Tranquada, C. D. Frost, 
R. J. Birgeneau, and H. Zhang, Phys. Rev. Lett. 98, 
247003 (2007). 

28 O.J. Lipscombe, S.M. Hayden, B. Vignolle, D.F. McMor- 
row, and T.G. Perring, Phys. Rev. Lett. 99, 067002 (2007). 

29 S.H. Naqib, J.R. Cooper, J.L. Tallon, R.S. Islam, and R.A. 
Chakalov, Phys. Rev. B 71, 054502 (2005). 



